Membrane buckling induced by curved filaments 
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We present a novel buckling instability relevant to membrane budding in eukaryotic cells. In this 
mechanism, curved filaments bind to a lipid bilayer without changing its intrinsic curvature. As 
more and more filaments adsorb, newly added ones are more and more strained, which destabilizes 
the flat membrane. We perform a linear stability analysis of filament-dressed membranes and find 
that the buckling threshold is within reasonable in vivo parameter values. We account for the 
formation of long tubes previously observed in cells and in purified systems. We study strongly 
deformed dressed membranes and their bifurcation diagram numerically. Our mechanism could be 
validated by a simple experiment. 



Eukaryotic cells are highly compartmentalized, and 
many of their confining structures are made of lipid bi- 
layers. In order to maintain the exchanges essential for 
their proper functioning, cells thus need tools that mod- 
ulate the shape and topology of these membranes. Such 
tools may be proteins that self-assemble to form tubes in 
solution pQ and can impose this intrinsic tubular shape 
on membranes j2]. In some other physically interest- 
ing cases, the structure of the protein does not suggest 
an obvious tubulation mechanism, as for the Endosomal 
Sorting Complex Required for Transport III (ESCRT- 
III) [3J. This protein complex is implicated in the for- 
mation of multivesicular bodies [I] , HIV budding [5] and 
cytokinesis [6 a , three processes which involve deforma- 
tion of the membrane into a bud and/or severing off the 
resulting membrane protrusion from the inside. Deep- 
etch electron micrographs of COS-7 cells overexpressing 
hSnf-7, one of the constitutive proteins of ESCRT-III, re- 
veal circular arrays of curved hSnf-7 polymers under the 
plasma membrane [Fig. [TJa)] [7]. This is evidence of the 
strong affinity of these filaments for the membrane [5] and 
for each other [9], as well as of their intrinsic curvature. 
When an ATP-hydrolysis deficient mutant of VPS4 — an 
ATPase involved in the disassembly of ESCRT-III fila- 
ments OHO] — is present, long membrane-covered tubes 
of hSnf-7 filaments are observed [Fig.[ljc)]. Similar struc- 
tures appear in in vitro systems using purified proteins 
[TT] , This suggests that tubes always form in vivo, but 
that in the presence of normal VPS4 alone they are im- 
mediately cut off the membrane to form vesicles. In this 
Letter we propose that this flat-to-tubular transition is 
a general feature of systems where curved filaments with 
attractive interactions bind to a membrane, and study 
this physical effect akin to the buckling of a rod (Fig. [TJ . 

We consider an infinite, initially flat lipid bilayer 
parametrized by its radial coordinate r. A subdomain 
Ti < r < r e of this surface is bound to an array of fila- 
ments (Fig. [TJ. The dressed membrane is then put into 
contact with the cytoplasm, which acts as a reservoir 
of filaments. In the following, we consider only axisym- 



metric configurations [T!5] and assume that the dressed 
membrane is very thin. We write the free energy of the 
dressed membrane as: 
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The first term is the Hclfrich free energy of the mem- 
brane, with bending modulus n, local total curvature c 
and tension a |12j . The second term represents the at- 
tractive interactions between filaments, characterized by 
a line tension 7. We assume that the filaments are closely 
packed; their surface density is thus constant through- 
out the array. The last term represents the free energy 
of the filaments. They have a preferred curvature r$ , 
and due to the cylindrical symmetry their actual curva- 
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Figure 1: (color online) Illustration of the proposed buckling 
mechanism. Overbent filaments are represented in blue, un- 
derbent filaments in red and the membrane in yellow. Wedges 
of the membrane were removed for visualization, (a) Curved 
filaments with an affinity for each other and the membrane 
form membrane-bound circular arrays. The tension and bend- 
ing modulus of the membrane tend to stabilize flat arrays, 
(b) Buckling, on the other hand, allows the binding of more 
filaments and the relaxation of those already bound to their 
preferred (yellow) radius. These stabilizing and destabilizing 
effects balance at the buckling threshold, (c) The formation 
of long tubes allows the binding of an arbitrarily large number 
of filaments close to their preferred radius. 
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ture is r _1 . A Taylor expansion about Tq 1 to second 
order yields the filament stiffness k. We denote by fi the 
chemical potential difference between hSnf-7 in the cyto- 
plasm and bound to the membrane. In vivo, the circular 
filaments pictured in Fig. [T] are not actually continuous 
and can be made of several consecutive shorter filaments, 
ft is also possible that the hSnf-7 filaments are "living" 
polymers and exchange monomers with the cytoplasm. 
Therefore we consider that filaments of any length are 
always available, and that their chemical potential per 
monomer does not depend on their length, so that fi is 
uniform throughout the filament array. We ignore holes 
in the array resulting from thermal fluctuation, which is 
correct in the limit of large binding energies HD]. Eq. 
also ignores the up-down asymmetry of the dressed mem- 
brane, a simplification discussed later. Finally, we define 
the scaled filament stiffness K — kj k and membrane ten- 
sion £ = (a — i^Tq/k + fe/(2/t). 

We first consider the stability of flat arrays of fila- 
ments [Fig. |l[a)]. In Ref. [7], it is observed that these 
arrays have a finite, rather well-defined external radius 
r e . We attribute this feature to a chemical equilibrium 
between hSnf-7 in the array and in solution. Minimizing 
T with respect to r e for a flat membrane (c = and 
<1A = 2irrdr), one finds that the array has a finite ex- 
ternal equilibrium radius only if k/2r 2 j > /x, i.e. only if 
it is more favorable for a filament to be in solution than 
bound to the rim of a very large (r e — > +oo) array. Line 
tension will shrink the array and make r e vanish unless 
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Under these assumptions we always have < r, < r e and 
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We now discuss the buckling of filament-dressed mem- 
branes [Fig. [TJb)]. Experimentally, it is observed that 
the typical length scale of a hSnf-7 protrusion is much 
larger than n and smaller than r e . We therefore assume 
for simplicity that = and r e — > +oo. We parametrize 
the dressed membrane by its altitude z(r). The equilib- 
rium states are the solutions of the force balance equa- 
tion jf^y = with boundary conditions ^(0) = and 

2^(+oo) = 0. Therefore, z(r) is defined up to an ar- 
bitrary additive constant. As in the case of a buckling 
rod [13], the buckling threshold is the set of parameters 
where non-zero solutions of the linearized force balance 
equation satisfying the boundary conditions exist. This 
equation reads: 
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where the primes denote differentiation with respect to 
the scaled radius R — r/u, and 
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The general solution of Eq. Q reads z'(R) = c±fi(R) 
c 2f2(R), where c\ and C2 are arbitrary constants and 



fl (R) = e - R ' 2 R^ K ' 2 U (-V, 1 + V4 + 2X, R) , 



f 2 (R) = e- R l 2 R^ 1+K ' 2 M (-V, 1 + V4 + 2K, R) (6) 



The confluent hypergeometric functions of the second 
kind U and M are defined in Ref. 14J. Non-zero so- 
lutions of this form satisfying the boundary conditions 
only exist for certain values of the parameters, thereby 
defining the buckling threshold. Two parameter regimes 
must be distinguished: 

- For v ^ N, we have the following asymptotic behaviors: 
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Thus /i diverges as R — > while Eq. ^ implies f 2 (0) = 
0. Hence the boundary condition z'(0) = imposes c\ — 
0. Similarly, ]i{R) diverges as R — > +oo, thus z'(+oo) = 
yields C2 — 0. Therefore there is no non-zero solution 
to the linearized buckling problem. 

- For v — n G N, the singular terms of Eq. ^ vanish 
and /i and ji are both proportional to the generalized 
Laguerre polynomials L^"\R) [T3]. Hence Eq. Q has a 
unique solution, up to an arbitrary amplitude C: 
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Since L„ (R) is a polynomial of degree n in R, z' n satisfies 
the boundary conditions for any n. Therefore, there is 
an infinity of buckling thresholds, one per integer v = n. 
This is again reminiscent of the buckling rod problem, as 
each normal mode z n of the dressed membrane has its 
own instability threshold (Fig. [2|. In the following we 
only consider the most unstable mode n = 0. 

We now study strongly deformed dressed membranes 
[Fig. [ljc)]. We first check that our model accounts for 
the existence of long dressed membrane tubes similar to 
those observed in Refs. [TJEBHIj. For a cylindrical pro- 
trusion of radius r* and length I 3> r tl one can neglect 
the rounded tip and base of the tube. Introducing a fic- 
titious vertical point force / pulling the membrane up at 
r = 0, we minimize the free energy Q = T — j '£ with 
respect to r t and £ and find 
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Figure 2: (color online) Normal modes z n (R) of the dressed 
membrane, (a) Spatial structure of the first four normal 
modes at their respective buckling thresholds for K — 2.5. 
(b) Thin black lines: buckling thresholds as a function of n, 
K and E. Thick black line: stability limit of long, cylin- 
drical dressed membrane tubes. Protrusions are obviously 
more stable at small E, where the destabilizing influence of 
the filaments overrides the stabilizing effect of the membrane. 
Therefore the nth normal mode of the flat dressed membrane 
is linearly unstable for parameter regimes located under the 
nth thin black line and long tubes exist only under the thick 
black line. Thin cyan (grey) line: parameter regimes com- 
patible with the experimental data of Ref. [7J. Symbols are 
referred to in the main text. 



Consider an equilibrium situation in which a long dressed 
membrane tube is held at a constant length by a force 
/ = ft- The force is then suddenly set to zero. In the 
case of an upward initial force ft > 0, the tube tends to 
retract. If f t < 0, on the contrary, I increases and the 
dressed membrane spontaneously tubulates. This corre- 
sponds to the region of Fig. [2jb) located under the thick 
black line. Interestingly, long tubes are always stable 
when the flat dressed membrane is linearly unstable, but 
the reverse is not true. Thus there exists a regime, lo- 
cated between the thick black line and the n = line of 
Fig.j2jTb), where the flat dressed membrane is metastable. 
This regime is compatible with biologically reasonable 
parameter values. Indeed, combining Eqs. ^ and Q9J) , 
one finds: 
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Inserting r t ~ 70 nm and r e ~ 200 nm [7J and the es- 
timates k = 20k B T, a = 10~ 5 N.m" 1 and 7 = 1 pN in 
this equation, we obtain a numerical relation between the 
scaled tension and filament stiffness characterizing the 
experiments of Ref. [7] . We plot this condition as a thin 
cyan (grey) line on Fig. [2jb). This line traverses both 
the metastable and unstable regions, making it possible 
that the experiments of Ref. [7] reflect either regime. 

We consider the possibility that the flat dressed mem- 
branes observed in Ref. |7j are indeed metastable. In this 
hypothesis, an important quantity is the energy barrier 
AT separating the flat state from the more stable, tubu- 
lated state. To compute AT, we numerically solve the 
full nonlinear shape equation of the tube on a finite do- 



main < R < 25 [TS]. Using Z = z/u, we define S as 
the arc length along the dressed membrane in the (R, Z) 
plane [Fig. |3ja)] . We parametrize the dressed membrane 
by R(S) and the angle ip(S) defined by 

R{S) = cosif>{S), Z(S) = -smip(S), (11) 

where the dots denote the differentiation with respect to 
S. Minimization of the free energy Q yields the shape 
equation of the dressed membrane: 
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where F — /u/(27tk). This equation is identical to 
Eq. Q in the small-?/; limit and to the bare membrane 
tube shape equation in the absence of protein (k = 0, 
fi = 0) [TS]. In the following, we discuss the specific ex- 
ample K = 2.5, but we believe that other values of K 
yield a similar behavior. Let us first comment on the 
three regimes presented in Fig. [3]ja-b). For S = 1.1 [in- 
dicated by A in Fig. |2jb)], tubes always retract in the 
absence of an external force, as shown in Fig.[3jb). Low- 
ering the surface tension to X = 0.89 (□), one reaches the 
boundary of the metastable region. For 0.39 < S < 0.89, 
a positive force is required to extract short tubes, but 
long tubes grow spontaneously unless opposed by a neg- 
ative F. At £ = 0.39 (O) an d lower, even short tubes 
grow spontaneously and can be maintained at a finite 
length only by a negative force. In Fig.[3^b), crossings of 
the horizontal axis by the force-extension curves denote 
solutions of the biologically relevant, F = problem, the 
stability of which is indicated by the sign of the curve's 
slope. Plotting the lengths of these protrusions as a func- 
tion of S, we obtain the diagram Fig. [3jc), where we ob- 
serve that the loss of stability of the n = mode studied 
above yields a subcritical bifurcation. Focusing on the 
metastable regime (0.39 < S < 0.89), we note that form- 
ing an infinitely long tube requires first extruding a short 
tube from the dressed membrane, which is energetically 
unfavorable. The associated energy barrier is given 
by the free energy of the unstable solutions represented 
by the main thin branch of Fig. [3jc). Integrating force- 
extension curves similar to those of Fig. [3jb), we calcu- 
late the work required to reach these solutions from the 
metastable, flat state and plot the results on Fig. [3^d). 
Under the effect of thermal fluctuations, an energy bar- 
rier of height AT is crossed at a rate T - 1 e - A - ?r /' c B T ; 
where r ~ ns is the characteristic relaxation time scale 
of the system. When AT is of the order of a few fc^T, 
thermal fluctuations are sufficient to ensure the buckling 
of the dressed membrane within experimentally observ- 
able time scales. This is, however, not the case here, 
and the large energy barrier makes thermally activated 
ESCRT-III-mediated budding extremely unlikely in most 
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Figure 3: (color online) Numerically computed mechani- 
cal properties of strongly deformed dressed membranes for 
K = 2.5. (a) Parametrization and profiles. (b) Force- 
extension curves (L = £/u). (c) Black lines: bifurcation di- 
agram for the F — problem. Cyan (grey) lines: changes 
induced by a weak asymmetry of the dressed membrane. In 
both cases, thick (thin) lines represent stable (unstable) so- 
lutions, (d) Activation energy a flat dressed membrane 
needs to reach the i — +oo buckled solution. 



of the metastable regime. Therefore, in vivo, ESCRT- 
III-mediated budding either takes place only in (or close 
to) the regime where the flat dressed membrane is lin- 
early unstable, or is assisted by some unknown active 
process (e.g. actin polymerization, which is regulated by 
the ESCRT-associated protein Alix [16 ). 

We now comment on two approximations used 
throughout this work. First, we assumed that the in- 
teractions between filaments and between filament and 
membrane are independent of the slope of the dressed 
membrane (i.e. of whether the filaments lie in the 
same plane or are stacked upon another). For small 
slopes, this dependence can be expanded as f-t(Vz) — 
Mo + ^2(Vz) 2 /2 + O [(Vz) 4 ] and yields the same lin- 
ear stability analysis as above provided we redefine £ = 
(<T — HQ — fi2) r o/ K +k/(2K). Second, we ignored in Eq. 
any terms violating the z — > — z spatial symmetry. These 
terms are allowed in general since the dressed membrane 
is not up-down symmetric, and might be responsible for 
the fact that buckling systematically occurs toward the 
outside of the cytoplasm [7] . Formally, such an asymme- 
try destroys the bifurcation studied here. If it is weak, 
however, a stable, almost flat configuration still exists for 
high tensions and loses stability close to the predicted 
£ = 0.39 threshold, as illustrated in Fig. |3jc). 

Finally, we believe that a better understanding of 
ESCRT-III-mediated budding could be gained by study- 



ing it experimentally in the absence of any active pro- 
cess. We propose an in vitro setup where an aspiration 
pipette is used to control the tension a of a giant unil- 
amellar vesicle [T7]- Introducing ESCRT-III proteins in 
the surrounding solution at a known concentration (and 
therefore at known /S) one could vary £ through a 
and directly measure the buckling threshold and its de- 
pendence on [i. Due to the existence of the metastable 
region, we also predict a hysteretic behavior. 

In this Letter we presented a novel buckling mechanism 
relevant for a wide range of systems involving interacting 
membranes and curved filaments (possibly including e.g. 
the one studied in Ref. [18]). Our robust qualitative and 
quantitative predictions can be tested in rather simple in 
vitro experiments and could shed light on the biological 
problem of ESCRT-III-mediated budding. 
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